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The effective potential $ of a classical ion in a weakly correlated quantum plasma in thermody¬ 
namic equilibrium at finite temperature is well described by the RPA screened Coulomb potential. 
Additionally, collision effects can be included via a relaxation time ansatz (Mermin dielectric func¬ 
tion). These potentials are used to study the quality of various statically screened potentials that 
were recently proposed by Shukla and Eliasson (SE) [Phys. Rev. Lett. 108, 165007 (2012)], Akbari- 
Moghanjoughi (AM) [Phys. Plasmas 22 , 022103 (2015)] and Stanton and Murillo (SM) [Phys. Rev. 

E 91, 033104 (2015)] starting from quantum hydrodynamic theory (QHD). Our analysis reveals 
that the SE potential is qualitatively different from the full potential, whereas the SM potential (at 
any temperature) and the AM potential (at zero temperature) are significantly more accurate. This 
confirms the correctness of the recently derived [Michta et ah, Contrib. Plasma Phys. 55, (2015)] 
pre-factor 1/9 in front of the Bohm term of QHD for fermions. 

PACS numbers: 52.65.-y, 52.25.Dg, 52.27.Gr 


I. INTRODUCTION 


Dense plasmas have recently gained growing interest 
due to their relevance for the interior of giant planets 
as well as for laser interaction with matter and inertial 
confinement fusion scenarios. Examples of recent exper¬ 
imental studies include the ultrafast thermalization of 
laser plasmas [1] or free electron laser excited plasmas [2] , 
inertial confinement fusion experiments at the National 
Ignition Facility [3] and magnetized Z-pinch experiments 
at Sandia [4]. Questions of fundamental theoretical im¬ 
portance are the conductivity and heat conduction, the 
energy loss of energetic particles (stopping power) in such 
a plasma, e.g. [5] or the temperature equilibration of the 
electronic and ionic components [2]. 

Despite recent advances in modeling and computer 
simulations a fully selfconsistent treatment of these, in 
general, highly nonequilibrium electron-ion plasmas has 
not been possible so far due to the requirement of the 
simultaneous account of electronic quantum and spin ef¬ 
fects together with the (possibly) strong ionic correla¬ 
tions. The main problem here are the vastly different 
time scales of electrons and ions resulting from their dif¬ 
ferent masses. A possible solution of this dilemma is a 
multi-scale approach that has been proposed by Ludwig 
et al. in Ref. [6] . It takes advantage of the weak electron- 
ion coupling that allows for a linear response treatment 
of the electrons. This idea has been used by Graziani 
et al. to decouple the electron kinetic equation using an 
STLS (Singwi-Tosi-Land-Sjolander) scheme [7] or a re¬ 
cently derived extension [8]. 

The key of this multiscale approach is to absorb the 
fast electron kinetics into an effective screened potential 

of the heavy ions with charge Q where the screening is 
provided by the electrons via a proper dielectric function 


e, e.g. m 
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taken in the static limit. We note that also the effect 
of streaming electrons maybe important in warm dense 
matter. Then the frequency argument becomes uj = ku 
which leads to wake effects, that are well investigated 
theoretically and experimentally for dusty plasmas, e.g. 
i9Hn] and quantum plasmas, e.g. Ref. [12] and references 
therein. However, this is beyond the scope of this paper. 

A similar but even simpler approach that is based on a 
quantum hydrodynamic model (QHD) and replaces the 
linear response potential of the ions by a simplified ex¬ 
pression that is derived from linearized QHD (LQHD), 
> <|)LQHD^ rppg expression for ^lqhd was de¬ 
rived by Shukla and Eliasson (SE) who predicted an at¬ 
tractive interaction between ions, even in the absence of 
streaming, Ue = 0 US]. Comparisons with density func¬ 
tional theory revealed that this is incorrect [T4|. This 
also underlined the limitations of QHD models to weakly 
correlated dense plasmas nann]. A recent overview and 
more references can be found in Ref. El. 

Recently, modified expressions for $lqhd derived 
by Akbari-Moghanjoughi (AM) [18] and Stanton and 
Murillo (SM) [19]. The latter is particularly interest¬ 
ing because it is also applicable to finite temperatures, in 
contrast to the SE and AM potentials that neglect ther¬ 
mal effects. The question arises whether the AM and SM 
potentials are more accurate than the SE potential and 
whether they allow for an extension of the applicability 
limits of QHD. 

It is the goal of the present paper to perform such an 
analysis. To quantify the accuracy of these three po¬ 
tentials we use, as a basis, the standard static poten¬ 
tials of weakly correlated electrons that can be derived 
from quantum kinetic theory [20] by linearization in the 
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FIG. 1. (Color online) Effective ion potential computed from 
the RPA dielectric function for four densities and T = 0 = 0. 
Inset shows an enlargement of the potential around its zero 
values to resolve the Friedel oscillations. 

external perturbation. This directly leads to the static 
potential 0 involving the random phase approximation 
(RPA) for the dielectric function and the Mermin dielec¬ 
tric function, respectively. Both include kinetic effects, 
finite temperature and, in the latter case, also collisions. 
Being rooted in a kinetic approach, the corresponding 
static potentials are, by construction, more accurate than 
the potentials ^lqhd thus, allow for a rigorous test 
of the validity of the latter. 

The main conclusions of our analysis are the follow¬ 
ing: 1.) the SE potential is qualitatively different from 
the kinetic results. 2.) the SM potential exhibits very 
good agreement with the RPA. 3.) the AM potential is 
in agreement with the SM potential, at T = 0, differ¬ 
ing only in the notation and 4.) the SM potential is in 
good agreement with the RPA result even at elevated 
temperature. This confirms the correctness of the re¬ 
cently reported [2T] correction factor 1/9 that has to be 
included in front of the Bohm term of the QHD equations 
for fermions. 

The paper is organized as follows: We first recall the 
RPA and Mermin dielectric functions in Sec. nn and dis¬ 
cuss their basic properties such as the reproduction of 
Friedel oscillations. After this we recall the SE, AM and 
SM potentials in Sec. ill and present numerical results 
that compare them to the potential computed from 
the RPA and Mermin dielectric functions. The paper 
concludes with a brief discussion of the physical origin 
of the factor 1/9. The main steps of its derivation are 
outlined in an appendix. 


II. EFFECTIVE ION POTENTIALS SCREENED 
BY THE RPA AND MERMIN DIELECTRIC 
FUNCTIONS 

The mean-field result for the electron dielectric func¬ 
tion is given by the random phase (RPA or Hartree or 
quantum Vlasov) approximation, enpAik^u), where cor¬ 
relation effects are neglected, and damping of collective 


oscillations is entirely due to Landau damping. Tem¬ 
perature effects are straightforwardly included by using 
the corresponding Fermi function. For an overview on 
the final expression, see the Appendix of Ref. [12] . The 
effective ion potential Q computed from the RPA di¬ 
electric function in the zero temperature limit is shown 
in Figs. 12 and for different densities. The zero temper¬ 
ature RPA potential reproduces the Friedel oscillations, 
as is illustrated in Fig. by comparison to the large r 
asymptotics, cos(2/ci7r)/r^, (red dashed line). While the 
oscillations are very shallow and hardly of relevance for 
the thermodynamics and transport properties in warm 
dense matter, the reproduction of the correct periodicity 
is a useful consistency and accuracy test of the Fourier 
transformation in Eq. 0 . 

For later reference we also provide the static long wave¬ 
length limit of the effective ion potential which is nothing 
but the Yukawa potential 0 

$y(r;n,T) = ^e-'=-^ (2) 

r 

with the familiar inverse Yukawa screening length, /cy, 
that interpolates between the Debye and Thomas-Fermi 
expressions, in the non-degenerate and zero temperature 
limits, respectively, 

/cy (n, T) =-/c|^^^^/^/_i/2(/3/i). (3) 

Here kpr = V^oOpIvp is the Thomas-Fermi wave num¬ 
ber, and 1-1/2 is the Fermi integral of order —1/2. Note 
that the potential 0 depends-via ky-on density and 
temperature. 

In the warm dense matter regime we expect the elec¬ 
trons to be weakly to moderately coupled so that corre¬ 
lation effects are, in general, not negligible. The simplest 
quantum dielectric function (DF) which takes collisions 
into account in a conserving fashion via a relaxation time 
approximation is the Mermin DF [22] 

eM{k,oj) = 1 + 

_ juj + iv) [eRPA {k,uj + w) - 1] _ 

uj + iu[eYiPAik,u) + iv) — l]/[eRPA(fc, 0) — 1] ’ 

which involves the finite temperature RPA DF and the 
electron collision frequency v. In Ref. m various choices 
for the collision frequency were studied, so here we will 
just use a typical set of values to highlight the main 
trends. 


III. TESTING STATICALLY SCREENED ION 
POTENTIALS DERIVED FROM QHD 

A. The Shukla-Eliasson (SE) potential 

The Shukla-Eliasson potential has been derived from 
linearized QHD, below we reproduce its most recent 
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FIG. 2. (Color online) Effective ion potential computed from 
the RPA dielectric function for two densities corresponding to 
Ts — 2.53 (top) and 0.25 (bottom) and T — 0 — {). The zoom 
into the small potential values confirms the correct oscillation 
frequency at large distances that is expected for Friedel os¬ 
cillations, as seen from the comparison to the (red) dashed 
curve. Here ro corresponds to the distance indicated by a 
vertical dashed line. The grey box shows that the oscillation 
period of the full potential agrees with the known asymp¬ 
totics. 

variant m- It involves the electron plasma frequency 
ujpe = (47rnoe^/em*)^/^, [m* is the effective electron 
mass] and the electron Fermi speed = fi(37r^)^/^/m*ro, 
where vq = ^ is proportional to the Wigner-Seitz ra¬ 

dius. The term = (0.328e^/m*ero)^/^ x [1+ 0.62/(1 + 
18.36a^ny^)]^/^ can be regarded as (optional) correction 
that is supposed to account for exchange-correlation ef¬ 
fects [it was proposed in Ref. [23] based on the exchange- 
correlation potential of density functional theory]. Since 
the other potentials do not include such a correction and 
it does not qualitatively alter the potential, in the com¬ 
parisons below we will use Vex = 0- 

The SE potential, for > 0.25, is given by [13]: 

= — [cos(/c^_r) b'sm{k'_r)] e~^+^, (5) 

with the parameter h' = 1/ The wave vectors 
are given by [24] k'^ = ks{V 4^*^^ ± l)^/^/\/4a^, where 
ks = ujpel^/vfJfT^ is similar (but not identical) to 
the inverse Thomas-Fermi screening length ky^ Eq. (©• 
The SE potential for the second parameter range. 


^SE ^ 0.25, is given by: 

(r; n) = [(1 + b) 6-'=+'- + (1 - 6) 6-'=-’’] . (6) 

2r 

with the definitions k± = ks{l T Vl — /V2a^^ 

and 6=1/ Vl — 4^*^^. 

Note that the Shukla-Eliasson potential is a zero- 
temperature approximation, and the potential depends 
only on density. Numerical results are included in Eig.[^ 

B. The Akbari-Moghanjoughi (AM) potential 

Here we reproduce the electrostatic potential derived 
by Akbari-Moghanjoughi (we use the most recent cor¬ 
rected version. Ref. HU). This potential has the same 
mathematical form as the SE potential , for > 
0.25: 

(/)^^>(r;n) = — [cos(/c'_r) + 6'sin(A:'_r)] e~^+^. (7) 

The difference lies in the definition of the parame¬ 
ter = Ji?ujp^l36m‘l{v‘^l3)‘^ and of the wave num¬ 

bers [24] k'^ = kTF{V4:a^^ ± l)^/^/V4(a^^ where 
krF = ^pe !is the inverse Thomas-Eermi screening 
length. These formulas are connected to the SE versions 
by '^ecc = 0 and an additional factor of 1/9 in the defini¬ 
tion of . 

The potential for < 0.25 is not given in Ref. m- 
Assuming the same relation to the SE result it has the 
form 

n) = ^ [(1 + b) e-'=+’’ + (1 - 6) , (8) 

with b = l/Vl — 4:a^^ and k± = T 

Vl -4(a^^)i/VV2(a^^. 

As the SE potential, this is a zero-temperature result 
which, therefore, depends only on the density. Numerical 
results are included in Eig. 

C. The Stanton-Murillo (SM) potential 

Here we reproduce the screened potential of Stanton 
and Murillo following the latest corrected version m- 
Again, there are two cases. Eor > 1 the potential is 
given by [25] 

{r;n^T) = — [cos(/c^_r) -f 6'sin(/c^_r)] e~^+'^^ (9) 

with = 3^/8PXILi/ 2 {rjo)/ it, A = 1/9, b' = 

— 1, and k'^ = kyEiVoi^^ ± /^/a^. 

Ip{r]) = f^dx x^/{l + denotes the Eermi inte¬ 

gral and Ip{r]) its derivative with respect to r]. r]o is de¬ 
termined by the normalization, no = \/2/i/2(^o)/7r^/5^/^ 
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with the inverse temperature /3. The inverse Thomas- 
Fermi screening length for finite temperatures is given as 
krp = (4/_i/2(%)A\/2^)^/^. 

A second version of the potential follows for the case 
< 1 : 

(r; n, T) = ^ [(1 + 6) 6-'=+’' + (1 - 6) e”''-’’] , 

( 10 ) 

where b = l/\/l — and k± = T 

Vl-a5M)l/2/y^W72. 

We note that the SM expressions are identical to the 
AM version in the limit T = 0, if is rescaled by a 
factor of 4, {T = 0). Numerical results are 

included in Figs. and 

D. Numerical test of the potentials $lqhd 
linearized QHD 

We now present a numerical comparison of the three 
LQHD potentials to the RPA potential 0- In the follow¬ 
ing, all results are given in Hartree atomic units. Since 
the AM potential conincides with the SM potential taken 
at T = 0 we will, in the following, not distinguish both 
potentials. We use the shortcut “SM” for both cases. 

1 . Zero temperature 

To begin the analysis, we recall that the SE, SM [and 
AM] ion potentials are all derived by applying the long- 
wavelength limit to the full potential. Correspondingly, 
in Fourier space one should expect accurate results only 
for small wave numbers, k 2k where kp = (Stt^u) is 
the Fermi wavenumber. 

For a comprehensive comparison we consider a broad 
density range, covering 6 orders of magnitude with 
the three values of the Brueckner parameter, rs = 
25.27,2.3,0.25 (similar values were studied in Ref. [E]). 
The first observation, cf. Fig. is that all four poten¬ 
tials show the same overall behavior with the deviations 
growing when the density is reduced. The next observa¬ 
tion is that the SM potential conincides with the RPA 
potential at large distances: for r/ap ^ 12,1.2,0.9, for 
the rs = 25.57, 2.3 and 0.25, respectively. However, a 
closer inspection of the large distance behavior shows 
that substantial deviations remain: none of the three ap¬ 
proximate potentials reproduces the Friedel oscillations 
of the RPA potential, see also Figs. and Neverthe¬ 
less, the SM potential provides a good fit to the RPA 
potential if one discards the (small) Fourier component 
with k = 2kp. 

Consider now the case = 2.3, cf. middle part of 
Fig-i Here one clearly recognizes the difference between 
the SM (and AM) potentials, compared to the SE poten¬ 
tial. While the SE potential exhibits strong deviations 
from the RPA potential, the SM potential is very close 




r/as 



FIG. 3. (Color online) Comparison of the Shukla-Eliasson 
(SE, blue dashes) and Stanton-Murillo (SM, red dots) [it co¬ 
incides to AM for T = 0] potentials to the zero-temperature 
full RPA potential (full black line) for three densities. Top: 
rs = 25.27, center: rs = 2.3, bottom: rs = 0.25. The SM 
potential shows increasing deviations from the RPA potential 
$ for decreasing density, cf. top and middle figures. At high 
density (bottom figure) the SM potential exhibits very good 
agreement with the RPA potential. The SE potential exhibits 
large deviations from the RPA at all densities. Insets show 
the behavior for large distances and confirm that all mod¬ 
els (except for the RPA) are unable to describe the Friedel 
oscillations correctly. 


to the RPA result. At the highest density, = 0.25, 
bottom figure, the deviations between SM and SE are 
very small. But this is not surprising as, in this case, the 
Bohm term in the QHD equations is negligible, so the in¬ 
fluence of the different prefactors of this term (1, for SE, 
versus 1/9, for SM) is insignificant. Einally, at the lowest 
density, = 25.27, top figure, all three potentials show 
substantial deviations from the RPA for small distances 
below lOa^, although the deviations of the SE potential 
are significantly larger, in particular for larger distances. 
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However, at these low densities, the plasma is strongly 
correlated and the RPA potential itself is not applicable. 

Thus, from the zero-temperature behavior we conclude 
that the SM (and AM) potential shows a substantial im¬ 
provement over the SE potential, confirming the correct- 
nesss and importance of the prefactor 1 /9 in front of the 
Bohm term [21]. In fact, we will show in Sec. IV that 
this coefficient is not a free parameter but, for the de¬ 
scription of the long-wavelength properties of the plasma 
(such as the statically screened potential) the value 1/9 
follows rigorously. 


2. Finite temperature 

Quantum hydrodynamics, as used in the plasma 
physics community, is a zero-temperature theory, so the 
results of SE [T3| for the effective ion potential and the 
recent correction by Akbari-Moghanjoughi [18] are re¬ 
stricted to the ground state. At the same time, for rele¬ 
vant applications to dense plasmas, warm dense matter 
or laser plasmas, finite temperature effects are usually 
non-negligible. It was shown in Refs. [niEi] that the 
QHD equations are directly linked to the Thomas-Eermi 
theory which allows for a straightforward incorporation 
of finite-temperature effects. This link was exploited in 
Ref. [19] to derive a statically screened potential for ar¬ 
bitrary temperatures. 

It is now interesting to compare the corresponding 
finite-temperature SM potential to the RPA potential at 
the same temperature. This is done in Eig. [^for differ¬ 
ent values of the degeneracy parameter 0 = kBT/Ep at 
constant density. The figure shows that, for high density. 
Vs = 0.5, the agreement is excellent, for all temperatures. 
This is not surprising since, at these densities gradient 
corrections (the Bohm term of QHD) are small, as noted 
above. Here the temperature dependence is governed by 
the one of the ideal Eermi gas (Eermi pressure). More 
interesting is the comparison at lower density, = 2.3, 
bottom part of Eig. where correlation effects (and the 
Bohm contribution) are significant. Here the agreement 
between SM and RPA potentials is slightly worse and 
improves with increasing temperature. This is due to 
smoothening of the Eriedel oscillations with increased 
temperature. 

Overall, we conclude that the agreement between the 
SM potential and the RPA, even at finite temperature, 
is very impressive. This confirms again that the SM po¬ 
tential correctly captures the long-wavelength properties 
of the static RPA potential 0 . 

Einally, let us analyze the effect of collisions (correla¬ 
tions) on the dielectric function. It is well known that 
these effects become increasingly important outside the 
weak-coupling limit, i.e. for > I. Yet it is not clear, 
a priori, how important these effects are for the static 
ion potential. We, therefore, computed the static poten¬ 
tial 0 with the Mermin dielectric function for a broad 
range of temperatures, densities and collision frequen- 
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FIG. 4. (Color online) Comparison of the Stanton-Murillo 
potential to the RPA potential for various temperatures at a 
fixed density. Top: = 0.5, bottom: = 2.3. As it is seen 

the SM potential accurately describes the screening effect at 
high temperatures and densities. 



FIG. 5. (Color online) Ion potential at a fixed density, 
Vs = 0.5 and two collision frequencies (in units of the plasma 
frequency). The behavior is representative for other densities 
and temperatures and indicates that the electron collisions do 
not affect the static screened potential, within the relaxation 
time approximation (Mermin model). 


cies. A typical result is shown in Fig. where the RPA 
potential is compared to the Mermin dielectric function 
with the collision frequency u = O.bujpi, at zero tempera¬ 
ture. The differences are extremely small, and the same 
behavior is observed at other densities and elevated tem¬ 
peratures. The explanation is, of course, that correlation 
effects are particularly relevant at short distances where 
the Coulomb potential is large anyway. We note that 
the situation is very different in streaming plasmas that 
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are out of equilibrium. Here scattering effects have an 
important impact on the effective ion potential [26] . 

IV. SUMMARY AND DISCUSSION 


(fermions) the same potential (i.e. 7 = 1) was postulated 
by Manfredi and Haas [36]. However, the correct value 
(at T = 0) is 7 = 1/9 which is clearly confirmed by the 
comparison of the SE and SM potentials to the RPA in 
the present work. 


A. Accuracy of the LQHD screened potentials 

We have presented an analysis of three recently pro¬ 
posed static ion potentials in quantum plasmas and com¬ 
pared them to the static limit of the RPA potential Q. 
The first of these potentials (SE) was the one derived 
by Shukla and Eliasson [13] who predicted an attractive 
minimum of the potential [cf. middle part of Eig.[^ that 
would be responsible for ordering of ions in dense plas¬ 
mas. Our analysis reveals that no such minimum exists in 
the RPA potential, so this minimum should be regarded 
an artefact of the used LQHD model m- 

The comparison with the two other potentials, the 
ones of Stanton and Murillo (SM) [19] and of Akbari- 
Moghanjoughi (AM) [18] which do not exhibit a compa¬ 
rable minimum reveals the origin of this discrepancy (see 
below). The SM (and AM) potential exhibits very good 
agreement with the RPA potential at T = 0 indicating 
that it correctly captures the long-wavelength properties 
of the RPA, in contrast to the SE potential. 


B. Implications for the Bohm potential for 
quantum plasmas 


Since the only difference between the SE potential (§ 
and the SM (and AM) potential § is the different pref¬ 
actor of the Bohm term in the underlying QHD equa¬ 
tions, the good accuracy of the latter indicates that the 
choice of the prefactor 1/9 is the correct one. In fact, this 
question is easily answered by recalling the one-to-one 
correspondence between QHD and Thomas-Eermi theory 
with gradient corrections [Hin]. In Thomas-Eermi the¬ 
ory and, more generally, in density functional theory, it 
is well known that an expansion of the (non-interacting) 
kinetic energy (together with exchange contributions) in 
a power series in terms of density gradients converges to 
the full Hartree-Eock energy [27]. In the present case, we 
start from the statically screened potential 0 in RPA, 
so we should expect agreement from a gradient expan¬ 
sion of the Thomas-Eermi kinetic energy T[n\ without 
exchange and correlation corrections. Since the Bohm 
term of QHD corresponds just to the first non-vanishing 
gradient correction, T 2 , (quadratic in density gradients, 
see Appendix), the pre-factors of both are linked insepa¬ 
rably via the equation [21] 


ST2[n-,'y] ^ ^ 

6n ' 8m 


Vn 

n 


- 2 - 




= 7Ve 


MH 


n], (11) 


where is the Bohm potential for quantum systems 

with Bose statistics. Eor the case of quantum plasmas 


C. Derivation of the first gradient correction to the 
kinetic energy for fermions. The coefficient 7 = 1/9 

Let us briefly discuss the origin and validity of the pref¬ 
actor 1/9. This, coefficient was first obtained by Kompa- 
neets [28] and Kirzhnitz [29], almost 60 years ago, who 
used a Matsubara Green functions method to obtain the 
gradient correction to the Thomas-Eermi energy func¬ 
tional. In the same year Golden m derived this factor 
by expanding the density matrix of a many-electron sys¬ 
tem. Later Hohenberg and Kohn [31] presented a sys¬ 
tematic derivation in terms of a gradient expansion to all 
orders, on the basis of the electron polarization function 
for zero temperature. Mermin [32] extended the Hohen¬ 
berg and Kohn method to finite temperature but did not 
give a derivation of the gradient correction at finite tem¬ 
perature. This was done by Krizhnitz et al. [33] and by 
Perrot [34] who minimized the grand canonical potential 
using the known form of the free energy of the electrons. 
At last, the Eeynman path-integral method was used by 
Yang [ 35 ] to obtain the single-particle Green function and 
to prove the factor 1/9. 

Over the recent decades occasionally additional deriva¬ 
tions have appeared. In particular, the recent analy¬ 
sis of Akbari-Moghanjoughi [18] essentially reproduces 
earlier derivations based on the calculation of the long- 
wavelength limit of the inverse of the Lindhard dielectric 
function, at T = 0, whereas the analysis of Stanton 
and Murillo [19] directly starts from Perrot’s result [34] . 

We also note that there exists a different approach to 
the gradient corrections in Thomas-Eermi theory where 
the prefactor of the respective term is used as a fit param¬ 
eter to reproduce improved simulations (that go beyond 
the RPA), e.g. [38]. Similarly, for the case of electrons 
bound in atoms, the prefactor can be optimized to bet¬ 
ter match the Hartree-Eock ground state energy, yielding 
7 = 0.2 [Ml[40]. However, in the context of the screened 
potential of free charged particles in a plasma one has to 
use the static long-wavelength limit of the inverse dielec¬ 
tric function which is exactly given by the RPA, so there 
exists no freedom of choice, see Sec. |rv^ 


D. Wavenumber expansion of the inverse 
polarization function and inverse dielectric function. 

Improved static ion potential 

A key input for the gradient correction of the kinetic 
energy is the wavenumber expansion of the inverse of the 
static density response function, i.e. of the longitudinal 
polarization function n(k, u; = 0). The result is (for the 
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derivation, see the Appendix) 
1 


i^(k) = 


2n(k) 


— do “1“ d2k + d4k + ••• 


( 12 ) 


This is a completely general result where the whole mi¬ 
crophysics is contained in the value of the expansion co¬ 
efficients. From this we can directly obtain systematic 
approximations for the statically screened ion potential 
0 without recourse to the connection between Thomas- 
Fermi theory and QHD, Eq. 0 - Indeed, what is needed 
to evaluate 0 is the inverse of the dielectric function 
which is related to the polarization function by [ 20 ] 


e(k,0) = l-^n(k). 


(13) 


From this and Eq. 0 we obtain the desired wavenumber 
expansion of the inverse static dielectric function [and of 
the screened potential ([T])] 


e-i(k,0) = 


‘2k^ (do T ...) 


(14) 


2 k‘^ (do + a 2 k‘^ + a/^k^ + ...)— 47 re^ 

Using, as the lowest order approximation, the RPA 
polarization, we obtain for the first five non-vanishing 
coeffcients (see Appendix) 


II 

0 

27re^ 

(15) 

a2{n,T) _ 
do{n,T) 

bi 

44 ’ 

(16) 

ai{n,T) 

&i - b2 

(17) 

ao(n,T) 

16kj, ’ 

de{n,T) 

-bf + 2&162 - bs 

(18) 

ao(n,T) 

64kl, 

dfi{n,T) _ 

bf — 3bfb2 + 62 + 26163 — 64 

(19) 

do{n,T) 

256A:|, 


where the bi involve Eermi integrals of different orders 
[see Eq. (43) of the Appendix)]. Note that the finite 


temperature results for d 4 ,d 6 and dg are given here for 
the first time. This result allows one to systematically 
derive approximations for the screened potential by con¬ 
sidering the long-wavelength limit, /c ^ 0 of the inverse 
dielectric function (14). 

Lowest order result: The lowest order is given by 
neglecting all /^-dependent terms in the parantheses of 
Eq. (14), d 2 = d 4 = ... 0, and we obtain 

= (20) 


k^+i 




Inserting this into Eq. 0 , immediately yields the 
Yukawa potential 0 . 

Second order result: The next order is obtained by 
retaining the terms quadratic in k, 


'(k,0) = 


t“(l+gfe") 


k‘^ + I 




+ ~^k^' 

do 


( 21 ) 



FIG. 6 . (Color online ) p refactor 7 in front of the Bohm po¬ 
tential cf. Eq. ( [11] ), for fermions as a function of tem¬ 

perature. In the zero-temperature limit, 7 approaches 1/9, 
and it increases monotonically to 1/3 [33]. In contrast, for 
bosons or spinless particles this factor equals 1 (dashed line), 
as in the case of a single quantum particle m- 


The result for d 2 is (see Appendix) 

ftV-3/2(^) 


02 {n,T) = 




( 22 ) 


is the approxima- 
and (for T = 0) by 


In the ground state, T = 0, this coefficient becomes 
a^\n,T) = — Equation ^1 ) 

tion used by Stanton and Murillo 
Akbari-Moghanjoughi (AM) [18]. The approximation of 
Shukla and Eliasson is obtained for T = 0 by using 
without the prefactor 7(0) = 1/9. This coefficient yields 
the ratio of the corresponding results for fermions and 
bosons (or spinless particles) [21] . Erom the temperature 
dependence of a 2 we can extract the temperature depen¬ 
dence of the coefficient 7 which is displayed in Eig. 

7 increases monotonically from 1/9 to 1/3 [33l |34] , for 
temperatures large compared to Ep- Evidently this fac¬ 
tor is crucial for the correct treatment of fermions within 
Thomas-Eermi theory or in quantum hydrodynamics. 

To compare the first and second approximations to the 
static inverse dielectric function, Eqs. (20) and ( [^ , we 
plot again the associated static potentials computed ac¬ 
cording to Eq. 0 . In the top two parts of Eig. we 
show again the zero temperature results of SE and SM 
and, in addition the Thomas-Eermi screended Yukawa re¬ 
sult, Eq. 0 . Thus, we are able to assess the quality of 
the different orders of the wavenumber expansion of 
in comparison to the full RPA result. Eor = 2.3 we see 
that the Yukawa potential (zeroth order of the expansion) 
provides a good approximation, but the next correction, 
proportional to that is contained in Eq. (21) and. 


thus, in the SM potential constitutes an improvement. 
In contrast, the SE potential which contains the k^ cor¬ 
rection with a nine times larger prefactor exhibits a poor 
performance and is even significantly less accurate than 

























the zeroth order correction alone. The same trend is ob¬ 
served at higher density (top figure) where the zeroth and 
first order are both indistinguishable from the full result 
whereas the SE potential is substantially less accurate. 

The bottom part of Fig.j^shows the performance of the 
zeroth and second order potentials [TF and SM, Eqs. (20) 
and ( [2T] ), respectively] at finite temperature. While at 
low temperature, 0 = 0.1, the second order provides a 
significant impovement, for higher temperature, 0 > 1, 
the expansion exhibits sign alternating convergence and 
both approximations are of comparable quality with a 
slightly better performance of the SM potential. In other 
words, for 0 > 1, the standard static Yukawa potential 
provides a fairly accurate description of the static limit 
of the RPA. 





FIG. 7. (Color onlii^ Top two figures: same as middle and 
bottom parts of Fig.j^but with the Yukawa potential, Eq. 
“TF”, added for comparison of the ground state behavior. 
Bottom: temperature dependence of the potentials. Same as 
bottom part of Fig. I^with the Yukawa potential § added. 

Higher order results. Correlations: It poses no 


principal problem to extend the above results for the in¬ 
verse dielectric function to higher powers of k. With the 
equations ([^ [T^ the next three higher orders are avail¬ 
able, for arbitrary temperature. This should improve the 
behavior of the static potential also at shorter distances. 
However, for the short-distance behavior also other ef¬ 
fects will be important, most noticeably exchanges and 
correlations. One way to include correlations is to use 
static local field corrections, as proposed in Ref. m- A 
more systematic and powerful approach is to use many- 
body theory which provides dynamic correlation correc¬ 
tions, e.g. Refs. m and m and the cited references. 


E. Conclusions 

The conclusion from the above comparison of approx¬ 
imate static ion potentials to the RPA result is that, in 
fact, there is no freedom in the choice of the prefactor of 
the Bohm term for quantum hydrodynamics for fermions. 
The QHD equations follow from the Thomas-Fermi the¬ 
ory by performing a systematic gradient expansion of the 
kinetic energy or, equivalently, a long-wavelength expan¬ 
sion of the inverse polarization function Eq. ( |l2[ ). 
This expansion directly leads to the coefficient 1/9, at 
T = 0, in front of the Bohm term compared to the cases 
of bosons or of a single-particle, whereas for higher tem¬ 
peratures the coefficient increases monotonically to 0.3 
(in ID and 2D different coefficients apply [21]). This co¬ 
efficient is a consequence of the Pauli principle and has 
important implications for a hydrodynamic modeling of 
dense quantum plasmas. For example, we have shown 
that a neglect of this coefficient gives rise to a statically 
screened ion potential that is far less accurate than the 
standard Thomas-Fermi (Yukawa) potential, rather than 
an improvement. It also indicates that earlier derivations 
of QHD equations are invalid for quantum plasmas and 
earlier QHD results that were using the prefactor 1 do 
not describe fermions, in general, nor dense plasmas, in 
particular. 
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APPENDIX: FINITE TEMPERATURE 
GRADIENT CORRECTIONS OF THE KINETIC 
ENERGY AND INVERSE POLARIZATION 


Here we present a systematic derivation of the gradient 
correction to the Thomas Fermi energy using Mermin’s 
approach for finite temperatures. This appears to be the 































most systematic and general approach and allows one to 
derive higher order corrections as well. 
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Gradient expansion of the kinetic energy 


Following Hohenberg and Kohn m, and Mermin [32] , 
the electron energy is a functional of the density: 

E[n\ = y*i;(r)n(r) dr+i drdF H-T[n], (23) 

where the first term is the energy due to external field 
and the second is the interaction energy of the electrons 
in mean field approximation. T[n] is the kinetic energy 
(we neglect additional exchange-correlation terms) which 
is well known in the local approximation (Thomas-Fermi 
theory), 

To[n] = J drt[n{r)]n{r), (24) 

where the kinetic energy density is related to the (local) 
Fermi energy, t[n] = To obtain corrections to 

this local approximation, we expand the kinetic energy 
in terms of the density perturbation, h = n(r) — no, 


T[n] = To [no] + JK{r — r')h(r)h(r') drdr'+ 

+ yT(r, r', r")h(r)n(r')n(r") drdr'dr" + .... (25) 

We now transform the second term to an expansion in 
powers of the density gradient (gradient corrections), by 
Fourier expanding the kernel {Q is the volume). 




(26) 


From Eqs. (25) and (26) we have, using the convolution 
theorem, 

T[n] = To[no] + ^ ^ K{k)h{kf + .... (27) 


We are now looking for the long-wavelength limit and 
expand 

^(k) = do + d2k^ H- d4k'^ + ..., (28) 

Since K{r) is real, Eq. (^) leads to ^(—k) = ^(k), 
therefore, the expansion (28) contains only even powers 
of k. Erom this, we obtain for K{r — r'), using ( [2^ , 

K{r - r') = (do - d 2 VV' + ...)d(r - r'). (29) 

Substituting ( j^ into ( [^ one obtains 


T[n] = To [no] + T 21 [n], (30) 

z=i 

T 2 i[n] = J dr a 2 i [n] \ V'n(r) p, (31) 


where we took into account that /d(r) dr = 0 and iden- 
tifiy a 2 = — 0 . 2 . 

Thus we have rewritten the kinetic energy in terms 
of a local term plus gradient corrections. It was shown 
in Ref. [ 21 ] that this energy functional leads directly, via 
functional differentiation with respect to the density pro¬ 
file, to quantum hydrodynamic equations with the lo¬ 
cal term yielding the Eermi energy and the first gradient 
term giving rise to a Bohm-type potential. 


ST2[n-,j] _ ^ 

Sn ^ Sm 


Vn 

n 


2 


2 ^ 

n 


= jV^^[n], (32) 


where is the potential postulated by Manfredi and 
Haas [36]. What remains is to establish the relation be¬ 
tween the coefficient 7 and the coefhent a 2 of the ex¬ 
pansion (31). In the following we prove that, for zero 
temperature, to lowest order in the correlation energy, 
a 2 = which means that 7 = 1/9. 


Finite temperature first order gradient term in RPA 


To evaluate the coefficient a 2 = — 0.2 we recall that the 
energy associated with a charged particle density fluctu¬ 
ation can be expressed, in linear response, in terms of an 
effective potential (external plus induced potential) 

T[n] = To[no] + T ^^(k)i/,ff(k). (33) 

k 


Recalling the definition of the longitudinal polarization 
function, n(k) = n(k)/l/eff(k), Eqs. (27) and (33) allow 
us to identify 


Kik) = 


1 

2n(k) ’ 


(34) 


and to make use of known results for H below. The low¬ 
est order many-body approximation for H is the random 
phase approximation (RPA) which reads, for an arbitrary 
temperature m 


7 2 2 

IIrpa {k,uj) = - [ 9 {u + z) - g{u - z)], (35) 

where u = uj/{kvF), z = k/{2kF)^ Xo — 3/16 {hwp/E f)‘^^ 
kF = (37r^n)^/^, ujp = j and 

s(.) = -,(-x) = 

(36) 

where 0 = ksT/EF^ and r] = ji/ksT is the chemical 
potential. To obtain the long-wavelength limit of A, 
Eq. ( [ 2 ^ , we now expand the inverse of the static (i.e. 
uj = u = 0) polarization function in the limit z 1: 


X Ey 
x-y 
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KnPAik) = - 


+ 


^H 2 {v) 


2k%xlHj{7^) 24kUlHUv) 


fc2, (37) 


Hiiv) = ^ I-iMv), H2 {v) = ^ I-3/2{v). 


(38) 


The first term in (37) yields the coefficient ho in Eq. 


which does not contribute to the kinetic energy since the 
integral of n vanishes whereas the second term in (37) 
yields the coefficient a 2 : 


ao (N ,0) = - 
Ci2 (N ,0) = - 


2 klxlHi{vy 

72menHf{r]) ‘ 


(39) 

(40) 


In the ground state, 0 ^ 0, we have the limits Hi{r]) 1 
and H 2 — —1, so H 2 IHI —1, and (39) gives a 2 ( 0 ) ^ 
h^/{12m) meaning that the coefficient in front of the 
Bohm potential in Eq. ( [3^ becomes 7 = 1/9. The tem¬ 
perature dependence of this coefficient was displayed in 
Eig.[§ 

Einally, substituting a 2 from (39) to the formula (pOl), 


for the kinetic energy at finite temperature we have 


TW=T.[..ol+i|^ 


Sm-e J I-i/2(v) n 


dr, 


(41) 


were we eliminated the Eermi integral of order —3/2 using 
the formula I- 3/2 = ~‘^^-ij 2 ' 


Higher order gradient terms in RPA for finite T 

Now consider the higher order (/ > 1 ) terms of the 
gradient expansion (30). In the long-wavelength limit, 


z < 1 the kernel K can be written as: 


7^rpa(^) = - 


^2/^20-1/2 


mefcFd-l/2(l + Ei biZ'^^) ’ 


(42) 


where 


bi{[n],e) = 


I-i-i/ 2 iv) 
2i + 1 I-i/2{'i]) 


(43) 


Expanding (42) it is straightforward to obtain the higher 
(even) order terms of the gradient correction: 


T 4 n] 


h^TT^ 

Idrrie 


Te[n] -- 

TsM = 


_ 7 ^ 

64me 

256me 


/ 

/ 


bl - &2 

0[n]^/^ 1 / 2 ( 7 ?) 


fef — 261&2 + bs 
0[n]^/^ kl,I_i/ 2 {ri) 


I V^n |2 dr, 


I V^n |2 dr, 


/■6t-362&2 + 6i + 2&i63-64 , „4 , 2 . 

/ -^ ^ dr, 

J 0[n]'/2fc9J_i/2(r?) 


In the zero-temperature limit, T 4 and Tq were obtained 
by Hodges im and Murphy [42], respectively. Eor finite 
temperature these terms and, additionally, Tg are given 
here for the first time. 


Here we do not discuss the convergence of this gradient 
series. Obviously, before using higher order terms one 
should verify that all terms of the same order arising 
from three-particle correlations [third term in (25)] are 
included as well. Eurther, it is well possible that there 
exist comparable terms arising from correlation effects 
beyond the RPA, e.g. [43] . 


Einally we note that these results apply only to three- 
dimensional plasmas. Results for ID and 2D systems 
were discussed in Ref. [ 21 ] . Extensions of these results to 
finite temperatures are given in Ref. [44]. 
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